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A SURVEY FOR CIRCUMSTELLAR DISKS AROUND YOUNG SUBSTELLAR OBJECTS 
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ABSTRACT 

We have completed the first systematic survey for disks around spectroscopically identified young 
brown dwarfs and very low mass stars. For a sample of 38 very cool objects in IC 348 and Taurus, we 
have obtained L'-band (3.8 fim) imaging with sufficient sensitivity to detect objects with and without 
disks. The sample should be free of selection biases for our purposes. Our targets span spectral types 
from M6 to M9.5, corresponding to masses of ^^15 to 100 Mj^p and ages of <5 Myr based on current 
models. None appear to be binaries at 0'.'4 resolution (55-120 AU). Using the objects' measured spectral 
types and extinctions, we find that most of our sample (77% ± 15%) possess intrinsic IR excesses, 
indicative of circum(sub)stellar disks. Because the excesses are modest, conventional analyses using only 
IR colors would have missed most of the sources with excesses. Such analyses inevitably underestimate 
the disk fraction and will be less reliable for young brown dwarfs than for T Tauri stars. The observed IR 
excesses are correlated with Ha emission, consistent with a common accretion disk origin. In the same 
star-forming regions, we find that disks around brown dwarfs and T Tauri stars are contemporaneous; 
assuming coevality, this demonstrates that the inner regions of substellar disks are at least as long-lived 
as stellar disks and evolve slowly for the first ^--^3 Myr. The disk frequency appears to be independent 
of mass. However, some objects in our sample, including the very coolest (lowest mass) ones, lack IR 
excesses and may be diskless. The observed excesses can be explained by disk reprocessing of starlight 
alone; the implied accretion rates are at least an order of magnitude below typical values for classical 
T Tauri stars. The observed distribution of IR excesses suggests inner disk holes with radii of >2Ri,, 
consistent with the idea that such holes arise from disk-magnetosphere interactions. Altogether, the 
frequency and properties of young circumstellar disks appear to be similar from the stellar regime down 
to the substellar and planetary-mass regime. This provides prima facie evidence of a common origin for 
most stars and brown dwarfs. 

Subject headings: stars: low-mass, brown dwarfs — infrared: stars 



1. INTRODUCTION 

Brown dwarfs are now being found in abundance. In 
the solar neighborhood, the advent of very wide-field sky 
surveys have led to the discovery of hundreds of nearby 
ultracool dwarfs, many of which are likely to be substellar 
(e.g. Kirkpatrick et al. 1999; Delfosse et al. 1999; Burgasser 
et al. 2002; Hawley et al. 2002). Discoveries of more dis- 
tant substellar objects are increasing as deep, blank sky 
surveys become more powerful (e.g. Cuby et al. 1999; Liu 
et al. 2002). Likewise, surveys of nearby open clusters, 
in particular the Pleiades, have been very successful at 
identifying brown dwarfs (e.g. Bouvier et al. 1998; Stauf- 
fer et al. 1998; Moraux et al. 2001). Searches of young 
(<10 Myr) star-forming regions have also been quite re- 
warding. Substellar objects are hotter and more luminous 
in their youth, and hence it is feasible to probe to very 
low masses in star-forming regions (e.g. Luhman et al. 
2000; Najita et al. 2000; Zapatero Osorio et al. 2000; Lu- 
cas et al. 2001; Liu 2001). Though mass determinations 
for such young objects are necessarily dependent on theo- 
retical models, the coolest objects seem to be <15 Afj„p, 



comparable in mass to planets around old solar-type stars 
found from radial velocity studies (e.g. Marcy & Butler 
1998). 

While the number of known brown dwarfs is growing 
rapidly, the origin of these objects is an unanswered ques- 
tion. The initial mass function (IMF) measured in young 
clusters seems to be continuous from the stellar to sub- 
stellar regime (e.g. Hillenbrand & Carpenter 2000; Najita 
et al. 2000), suggesting brown dwarfs form in the same 
manner as stars. However, the overlap in mass between 
young "planetary-mass" brown dwarfs (^15 Mj^p) and 
planets discovered by radial velocity surveys is intriguing, 
maybe suggestive of a common genesis in disks around 
young stars. However, it is unlikely such brown dwarfs 
could form like giant planets are believed to, namely core 
accretion followed by runaway growth (Pollack et al. 1996). 
Perhaps gravitational disk instabilities (Boss 1998; Pick- 
ett et al. 2000) or disk-disk collisions (Lin et al. 1998; 
Watkins et al. 1998a, b) are involved, followed by ejection 
into the field. Indeed, a few solar-type stars are known 
to have massive substellar companions (Els et al. 2001; 
Marcy et al. 2001; Liu et al. 2001; Potter et al. 2002; Udry 
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et al. 2002), which might point to concurrent formation 
of brown dwarfs and planets. An alternative formation 
scenario involving ejection has been proposed by Reipurth 
& Clarke (2001). They suggest that brown dwarfs form 
as stellar embryos in small newborn multiple systems, but 
their growth is prematurely truncated by ejection due to 
dynamical interactions at very young ages ('~10'' yr). A re- 
cent computer simulation of cloud fragmentation by Bate 
et al. (2002) portrays an amalgam of these scenarios: most 
brown dwarfs form via gravitational instabilities in circum- 
stellar disks and then are ejected before they can accrete 
to stellar masses. 

The abundance of isolated brown dwarfs in star-forming 
regions and open clusters stands in stark contrast to the 
dearth of 20-70 Mj^p objects found as close radial ve- 
locity companions to mature, solar-type stars, a.k.a. the 
"brown dwarf desert" (Campbell et al. 1988; Marcy & 
Benitz 1989; Walker et al. 1995; Halbwachs et al. 2000; 
Zucker & Mazeh 2001; Pourbaix & Arenou 2001). But 
the implication of this phenomenon is unclear. It might 
suggest that isolated brown dwarfs and massive planetary 
companions arise from independent formation pathways. 
Alternatively, the difference might stem from evolutionary 
effects which preferentially deplete brown dwarf compan- 
ions, such as inward migration into the central star (Ar- 
mitage & Bonnell 2002) or outward ejection. 

One potential insight into the formation mechanism(s) 
for brown dwarfs is whether young substellar objects pos- 
sess circumstellar disks. There is abundant observational 
evidence for accretion disks around young solar-type stars. 
Moreover, disks are a key element in current theoretical 
understanding of star and planet creation: the formation 
of a disk is believed to be an inevitable step, since the in- 
falling natal material has angular momentum. Thus, the 
presence of disks around young brown dwarfs is naturally 
accommodated in "star-like" formation scenarios. On the 
other hand, scenarios involving collision and/or ejection 
are likely to be hostile to circumstellar disks. Hence, de- 
termining the frequency and physical properties of disks 
around young brown dwarfs, as well as their relation to 
disks around young stars, is an important observational 
goal. 

Evidence for circumstellar disks around individual 
young (<few Myr) brown dwarfs has recently been found. 
Strong Ha emission, likely due to ongoing accretion, has 
been seen from sources in Taurus (Bricefio et al. 1998; 
Martin et al. 2001), IC 348 (Luhman 1999), p Oph (Luh- 
man et al. 1997), and a Ori (Barrado y Navascues et al. 
2001; Zapatero Osorio et al. 2002a, b). Direct kinematic 
evidence for accretion from asymmetric Ha emission line 
profiles has been found for one object in Taurus at the 
substellar mass limit (Muzerolle et al. 2000). Infrared (IR) 
excesses have been detected from JHK (1.1—2.4 p,m) data 
for a few spectroscopically identified young brown dwarfs 
in p Oph (Wilking et al. 1999; Gushing et al. 2000), IC 348 
(Luhman 1999), and the Trapezium cluster (Hillenbrand 
et al. 1998; Lucas et al. 2001). A recent JHK survey by 
Muench et al. (2001) finds a majority of photometrically 
selected brown dwarf candidates in Trapezium have ex- 
cesses. 

IR excesses are one of the classic signatures for disks 
around T Tauri stars. Such excesses are correlated with 



UV/optical excess emission, as expected for an accretion 
disk origin (Hartigan et al. 1990; Valenti et al. 1993; Ed- 
wards 1995; Kenyon & Hartmann 1995). In particular, 
there is a one-to-one correspondence between IR color ex- 
cesses and optical veiling indicative of accretion onto the 
central source (Hartigan et al. 1995). There is other plen- 
tiful evidence for disks around T Tauri stars, including far- 
IR photometry, millimeter wavelength studies, and high- 
resolution optical/IR imaging — for nearly all objects, the 
different disk indicators agree, which supports the accu- 
racy of using a single indicator, in this case IR excess. 
On physical grounds, disks are the best explanation for 
both the IR excesses and the relatively low extinctions to 
the central sources (e.g. Adams & Shu 1986; Myers et al. 
1987). This body of evidence natiu'ally explains the IR 
excesses detected from young brown dwarfs as originating 
from circumstellar disks. 

Mid-IR (7—14 /im) detections are available for young 
brown dwarfs in p Oph (Comeron et al. 1998; Bontemps 
et al. 2001; Testi et al. 2002; Natta et al. 2002) and 
Chamaeleon (Comeron et al. 2000; Natta & Testi 2001; 
Apai et al. 2002), which further support the disk hypoth- 
esis. The mid-IR fluxes arc >10x greater than expected 
from the photosphere, and the resulting spectral energy 
distributions (SEDs) are much broader than can be ex- 
plained by a non-disk origin. 

However, it is difficult to determine the frequency of 
disks around brown dwarfs from studies to date due to a 
combination of small number statistics, sample selection 
inhomogeneity, and, most importantly, choice of wave- 
length. A priori, brown dwarf disks are expected to be 
harder to detect than disks around stars because of lower 
contrast. Subst{41ar objects are less luminous and have 
shallower gravitational potentials; hence, the inner regions 
of their disks are likely to be cooler and thus could have 
negligible excesses in the JHK bands, which have been 
used by most previous studies. Disks are more luminous 
at longer wavelengths, and they can be readily identified 
at mid-IR wavelengths. But because such observations 
can only detect young brown dwarfs with large excesses, a 
complete census of the disk frequency is not yet possible 
in the mid-IR. 

L'-band (3.8 pva) observations are an excellent means 
to address this question. Disks can be readily identified 
from ground-based telescopes by excess L'-band emission, 
which arises from warm circumstellar material within a 
few stellar radii (< 0.1 AU). In addition, existing infrared 
cameras can detect young brown dwarf photospheres, and 
hence the absence of a disk can be discerned. Such obser- 
vations have been used to measure the disk frequency for 
low-mass stars in young clusters (e.g.. Haisch et al. 2000). 
These studies find a higher disk frequency than JHK or Ha 
studies, as expected from the greater sensitivity of thermal 
IR observations to less luminous disks. 

In this paper, we present results from a large L'-band 
survey for circumstellar disks around young brown dwarfs 
and very low mass (VLM) stars in the IC 348 and Taurus 
star forming regions. The incidence and properties of disks 
around young substellar objects may shed light on the ori- 
gin of these intriguing objects. Our study targets young 
sources which have been spectroscopically classified to be 
very cool, with spectral types of M6 and later. By focus- 



Disks around Young Brown Dwarfs 



3 



ing on targets with spectral types, we are more sensitive 
to small IR excesses since we can determine the intrinsic 
photospheric colors. Our selection criterion also ensures 
that we are targeting very low mass members, lying near 
or below the stellar/substellar mass boundary. Since our 
observations are sensitive enough to detect bare photo- 
spheres, this is the first unbiased survey for disks around 
such objects. 

We describe our sample and our observations in § 2. We 
discuss the frequency of IR excesses and correlations with 
other properties in § 3. Simple models are used to study 
the disk characteristics in § 4. We examine the possible 
systematic effects in § 5. We then discuss the implications 
of our findings for brown dwarf disk searches, the proper- 
ties of such disks, and brown dwarf formation scenarios in 
§ 6. Finally, we summarize our results in § 7. For readers 
primarily interested in the results, we suggest focusing on 
§ 6 and § 7. 

2. OBSERVATIONS 

2.1. Sample Selection and Properties 

For objects in young star-forming regions, the substellar 
boundary is believed be around spectral type M6, based on 
comparison with theoretical models (Luhman et al. 1998a) 
and the application of the lithium test in young clusters 
(Basri 1998). Most young brown dwarf candidates do 
not have lithium measurements, which would guarantee 
their status as cluster members. However, their positions 
in color-magnitude diagrams indicate they arc very likely 
cluster members and the very late spectral types are con- 
sistent with substellar status. 

Our sample of 38 objects includes nearly all published 
objects with spectral types of M6 or later in the IC 348 
and Taurus star-forming regions. We adopt a distance of 
140 pc to Taurus (Elias 1978; Kenyon et al. 1994) and 
300 pc to IC 348 (Herbig 1998, and references therein), 
with an error of 10% for each. The typical stellar ages in 
these regions have been estimated to be around 1—3 Myr 
(Kenyon & Hartmann 1995; Herbig 1998; Palla & Stabler 
2000). There are three components to our sample: 



IC 348 core region: these objects originate from the 
NICMOS imaging survey of Najita et al. (2000). 
By using narrow-band photometry to measure the 
depth of the 1.9 /im H2O absorption band, they de- 
rived spectral types and extinctions for all objects 
down to K = 16.5 mag in the central « 5'x5' re- 
gion.^ This sample of substellar objects is unique 
because it is complete, in the sense that every ob- 
ject in the survey region has an estimated spectral 
type. Among their identified cluster members, we 
observed nearly all objects with spectral types from 
M5.7 to M9.4, with the exception of two (objects 
064-05 and 014-03, with spectral type of M5.8 and 
M6.3, respectively). Cluster members were selected 
as objects being 5 Myr or younger in their color- 
magnitude diagram. We also observed one object 
(071-01, spectral typo M5.8) with an estimated age 



of wlO Myr. A number of these sources have pub- 
lished optical spectroscopy (see § 5.1), and we have 
also been pursuing ground-based optical and near- 
IR spectroscopy (M. Liu et al., in preparation); the 
spectroscopy confirms that these sources arc gen- 
uine very cool (late M-type) young objects. Since 
this NICMOS sample has been selected and ana- 
lyzed in a homogeneous fashion, our analysis will 
often examine the properties of this sub-sample in 
addition to those of our sample as a whole. 

• IC 348 outer region: we include additional very cool 

members of IC 348 from Luhman (1999), originally 
selected from R and /-band optical imaging. These 
sources lie outside the Najita et al. (2000) survey 
region. Luhman (1999) used far-red optical spec- 
troscopy to derive spectral types and J-band extinc- 
tions. No errors are given for his extinction determi- 
nations; we conservatively assiime these are 0.1 mag. 

• Taurus: we include nearly all published objects with 
spectral types of M6 or later in the Taurus star- 
forming region. These objects were originally found 
from optical or X-ray surveys. Four of them come 
from Martin et al. (2001), who derived spectral 
types from far-red optical spectroscopy and V-band 
extinctions from broad-band colors. The remain- 
ing objects have spectral types from Briceno et al. 
(1998) and Luhman et al. (1998a); for these, we 
use ff-band extinctions derived by Luhman (2000) 
from broad-band colors. No errors are given for the 
published extinction determinations; we conserva- 
tively assume these are 0.5 mag and 0.1 mag for the 
Martin et al. (2001) and Luhman (2000) results, re- 
spectively. 

Figure 1 plots our targets in comparison to theoretical 
models from Baraffe et al. (1998, 2002), D'Antona & Mazz- 
itelfi (1997), and Burrows et al. (1997). For the D'Antona 
& Mazzitelli models, we use a newer unpublished version, 
which we refer to as the DM98 models.^ To compare the 
models to the observations, we use a linear fit for the con- 
version between spectral type and effective temperature 
(Toft ) advocated by Luhman (1999) for pre-main sequence 
objects: 

T^fi{K) = 3850 - 141.0 x SpT (1) 

where SpT is the M spectral subclass, i.e., SpT = means 
MO, SpT ^ 5 means M5, etc. To plot the DM98 and Bur- 
rows et al. models, we also use a linear fit to if-band bolo- 
metric corrections derived by Leggett et al. (2000, 2001): 

BCk = 3.99 - Teff/2424 . (2) 

iO.l mag in BCk, negligible 



The scatter about the fit is 
for our purposes. 

The model-derived masses range from Ril5 to 100 Mj„p. 

These depend on the choice of models and temperature 
scale; our choices lead to conservative (i.e., larger) values 

^ The spectral types in the IC 348 core sample are estimated from the correlation between /-band spectral types and the strength of the 1.9 /um 
water band, as measured for both field dwarfs and a subsample of young stars in IC 348 itself (see Najita et al. 2000 for details). The use 
of fra<;tional M subtypes reflects the continuous nature of the water index measurement and is not meant to advocate a new spectral typing 
convention. 

^ The models are available at http: //www. mporzio. astro. it/"dantona. 
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for the masses (see § 5.5). The ages inferred from the mod- 
els arc <5 Myr, and do not depend strongly on the choice 
of temperature scale since the isochrones are roughly hor- 
izontal in the HR diagram. Taken at face value, some 
models indicate that the lower mass objects tend to be 
younger than the higher mass ones. However, given the 
substantial uncertainties in the models of young substel- 
lar objects (Baraffe et al. 2002), we largely avoid using 
model-derived quantities in our analysis. Instead, we as- 
sume that the stars and brown dwarfs in these regions are 
approximately coeval, and adopt the age estimate of a few 
Myr derived for the stars. 

2.2. L' -hand Imaging: UKIRT 

We obtained L'-band (3.4-4.1 /xm) imaging on 13—15 
October 2001 UT at the 3.8-m United Kingdom Infrared 
Telescope (UKIRT). We used the facility camera IR- 
CAM/TUFTI, which employs a 256 x 256 InSb array and 
has a pixel scale of O'.'OSl pixcl^^. The filter used was the 
L'-band filter produced by the Mauna Kea Filter Consor- 
tium (3.4—4.1 fjLm; Simons & Tokunaga 2001; Tokunaga 
et al. 2001). UKIRT has a fast-steering secondary mirror 
which provides tip-tilt correction by sensing a star off-axis 
of the science target. The typical resulting image quality 
was excellent (0'.'35— 0'.'4 FWHM), and none of our sources 
were resolved into binaries. Conditions were photometric. 
During the course of the night, we regularly observed stan- 
dard stars from a UKIRT in- house calibrator list. 

We used a standard "two-beam" observing technique 
in order to accurately cancel the very large and variable 
thermal flux from the atmosphere. For each target, we ob- 
tained a coaddcd set of short integrations, then offset the 
telescope by about half of the array and integrated again. 
This was repeated, with small dithers at each nod posi- 
tion to ensure the object did not fall on the same position 
on the array. Integration times were determined using the 
published K-h&Tid magnitudes, spectral types, and extinc- 
tions to reach comparable S/N for all the objects. The 
on-source integration ranged from 3 to 60 min. 

Reductions were done in standard fashion for thermal 
IR imaging. For each object, a flat-field was created from 
the average of all the images. Each flat-fielded image was 
then subtracted from the corresponding image in the other 
nod position. The images were registered using the relative 
telescope coordinates and averaged together. The result- 
ing mosaic contained a positive and negative image of the 
target. Both were measured, and the standard error of the 
results was adopted as the uncertainty. For the standard 
stars, we used a 5" radius, as is done for the UKIRT stan- 
dard star programme. We then determined a zeropoint 
and extinction coefficient for each night. Since most of 
our targets are faint, we used a smaller aperture (10 pixel 
radius) for their photometry, and then applied an aperture 
correction out to 5" based on curves of growth measured 
nightly from the much brighter standard stars. The size 
of this correction was on the order of 0.15 mags. 

Table 1 presents our resulting measurements. The final 
errors for the object magnitudes are the quadrature sum 
of the errors in the raw photometry (the dominant term), 
the zeropoint, the extinction coefficient, the aperture cor- 
rection, and the standard star magnitudes. 

One object in the photometry sample, 045-02 in IC 348, 



was significantly bluer in Ks — L' than the colors of late 
M dwarfs. Based on JHKL' colors, it is probably a back- 
ground late-type giant whose position in the observed 
color-magnitude diagram overlaps the cluster members. 
We include its photometry in the Table 1, but the object 
is excluded from the analysis. 

2.3. IR Photometry: 2MASS 

We applied very small transformations to all the rele- 
vant photometry, namely the IR colors for field and young 
cluster M dwarfs, to bring them to a common photomet- 
ric system. We chose to use the well-understood 2MASS 
system (Carpenter 2001). For the JHK colors of M dwarfs 
in the solar neighborhood, we used measurements from 
Leggett (1992, hereafter L92), Leggett et al. (1998, here- 
after L98), and Leggett et al. (2002, hereafter L02), ex- 
cluding the subdwarfs. The L92 and L98 near-IR data 
are on the CIT system, which we converted to the 2MASS 
system using Carpenter (2001). The data from L02 are 
on the MKO-NIR system, which wc converted rising addi- 
tional transformations from Hawarden et al. (2001). The 
transformed samples showed no significant differences in 
their near-IR colors. For the L'-band data, all the data 
were either in the UKIRT system (L92 and L98 data) or 
the MKO-NIR system (L02 data and our observations). 
These are essentially identical for M dwarfs so no trans- 
formations were needed (L02). 

The Najita et al. (2000) sample has been observed in 
the HST/NICMOS F166N and F215N filters, which arc 
1% bandwidth. These filters arc placed in feature-free 
regions of H- and K-hand stellar spectra. To convert 
these measurements to a standard photometric system, we 
matched 77 objects in the Najita et al. sample with objects 
in the 2MASS Second Incremental Data Release (Cutri 
et al. 2000). We found (F166N-i?) = 0.00±0.02 mags and 
(F215N-i4's) = 0.08±0.01 mags with no statistically sig- 
nificant dependence on J ~ K color. We added J-band 
data from 2MASS, available for 23 out of 25 objects. 

For the Luhman (1999) IC 348 sample, wc used 2MASS 
near-IR photometry, available for all the objects except 
one (#405). For that one, we used his published photom- 
etry. The published near-IR photometry for the Taurus 
sample all came from 2MASS. 

3. OBSERVATIONAL RESULTS 

We identify IR excesses around our targets using Kg — L' 
colors, which track the shape of the SEDs from 2.0-4.1 /zm. 

This choice is advantageous for several reasons, includ- 
ing: (1) disks arc more luminous at longer wavelengths, 
as compared to just using JHK data; (2) the effects of 
extinction (and uncertainties therefrom) are reduced at 
longer wavelengths; (3) the small wavelength range min- 
imizes the effect of extinction errors; (4) previous stud- 
ies of disks around T Tauri stars have used these wave- 
lengths; (5) color measurements are distance independent; 
and (6) there is ample published Kg — L' data for field 
M dwarfs. Note that all the results discussed in this sec- 
tion are strictly empirical^ e.g., independent of theoretical 
models and choice of temperature scale. 

3.1. IR Excess Frequency and Amplitude 
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The top panel of Figure 2 shows our observations: 
Ks — L' colors for the sample as a function of spectral type 
compared with the colors of the field M dwarf sample de- 
scribed in § 2.3. A substantial number of objects show 
Kg — L' colors redder than expected from purely photo- 
spheric emission, while the lower envelope of the observed 
color distribution agrees well with the dwarf locus. 

Since all objects in our sample have measured spectral 
types and extinctions, we can determine their intrinsic 
photospheric Kg — L' colors, and hence we are very sen- 
sitive to the presence of excess IR emission, i.e., the pres- 
ence of circumstellar disks. The bottom panel of Figure 2 
shows the dereddened Kg — L' colors, using the published 
extinctions and the extinction law from Mathis (2000). 
A few trends are apparent: (1) most objects have dered- 
dened Ks — L' colors in excess of that expected from their 
photospheres alone; (2) the lower envelope of the color 
distribution is consistent with the locus of field M dwarfs, 
suggesting that such objects provide a legitimate compar- 
ison (see also § 5.3); and (3) the maximum amplitudes of 
the IR excesses decrease at later spectral types, i.e., lower 
masses. 

For each object, we compute the intrinsic Ks — L' excess 
in the usual fashion: 

E{Ks-L')o = {Ks-LX"""""^ - {Ks-L%'""°'P^^"^ . 

(3) 

The photospheric colors as a function of spectral type are 
determined from the field M dwarf sample. To compute 
the error in E{Ks~' L')o, we add in quadrature the errors 
in the Ks — L' colors, the reddening, and the spectral types 
(which add an uncertainty to the intrinsic photospheric 
color). The resulting excesses and their errors are listed 
in Table 1. In general, the excesses are much smaller than 
those observed for T Tauri stars (e.g. Kenyon & Hartmann 
1995). 

Table 2 gives the breakdown of the IR excess fraction 
by spectral type and target subsample. Overall, 31 out of 
38 (82%) of the sample show {Ks — L')o excesses, with 21 
out of 38 (55%) having excesses larger than their la mea- 
surement error. To compute the disk frequency and its un- 
certainties, we use a Monte Carlo technique described in 
Appendix A which accounts for both the (gaussian) mea- 
surement errors and the (Poisson) counting errors. We 
find an excess frequency of 77% ± 15% for the sample — 
disks around young brown dwarfs and VLM stars appear 
to be very common. 

Next, we examine any correlations between the IR ex- 
cesses and other physical properties of the objects. Note 
that since the detected amount of disk emission depends 
on the viewing angle to the observer, any such trends are 
not expected to be very strong. Random variations in 
the viewing orientation will inevitably act to obscure such 
trends. 

3.2. Correlation with Ha 

For T Tauri stars, both the optical line emission and the 
IR excess are believed to originate from an accretion disk. 
The IR excess comes from warm dust grains in the disk, 
while Ha emission originates from the accretion of disk ma- 
terial onto the central star, e.g., via a boundary layer or a 
magnetically regulated accretion flow. Therefore, the Ha 
emission and IR excesses should be correlated, and indeed 



such a correlation is seen among T Tauri stars (e.g. Valenti 
et al. 1993; Kenyon & Hartmann 1995). Figure 3 shows 
such a comparison for our sample of young VLM stars and 
brown dwarfs, using Ha data mostly from the literature 
(Briceno et al. 1998; Luhman et al. 1998a, b; Luhman 1999; 
Martin et al. 2001; K. Luhman 2002, priv. comm.). The 
entire Taurus sample has published Ha measurements, and 
about half of the IC 348 sample does. A correlation is 
seen between E{Ks — L')q and Ha equivalent width, with 
a Spearman rank correlation coefficient rg = 0.64. The 
probability of this being drawn from a random sample is 
0.0009, or a 3cr correlation. This is comparable to the cor- 
relation observed for T Tauri stars (Kenyon & Hartmann 
1995) and provides strong circumstantial evidence for ac- 
cretion disks around young brown dwarfs. 

3.3. N on- Correlations with Mass and Age 

Figure 4 presents a histogram of {Ks — L')q excesses as 
a function of spectral type. For the three earliest spectral 
type bins (M5.7-M6.4, M6.5-M7.4, and M7.5-M8.4), a 
Kolmogorov-Smirnov (K-S) test (Press et al. 1992) finds 
a high probability (>40%) that all the samples are drawn 
from the same parent population. Since the evolutionary 
tracks at constant mass are roughly constant in Tcs for 
young ages, the spectral types correspond to a relative 
mass scale. Hence, we find no evidence for a strong de- 
pendence of IR excesses on mass. The exceptions are the 
coolest (lowest mass) objects, types M8.5— M9.4, where 
the excesses are small, consistent with non-existent. The 
disk frequency likewise appears to be independent of mass 
(Table 2), though larger samples of the coolest objects are 
needed to better study this issue. Using only the homoge- 
neous IC 348 core sample gives comparable results. 

In Figure 4, there is a hint that the Taurus population 
might tend to have larger {Ks — L')o excesses than the 
IC 348 samples. However, a K-S test of the M6-M7 ob- 
jects shows the Taurus {Ks — L')o excess distribution dif- 
fers from that of the IC 348 core sample at only the 1.5a 
level. Also it is not clear if/how the selection eff'ects of the 
Taurus sample impact such a comparison. Hence, we find 
no strong evidence for differences in the sub-samples' IR 
excesses. 

Figure 5 examines the age dependence of the inner disk 
properties by plotting the {Ks — L')q excesses as a func- 
tion of dereddened absolute X^-band magnitude, Mks- 
If we separate objects by spectral subclass, Mks can be 
used as a surrogate for age, due to the fact that model 
isochrones are roughly horizontal in the HR diagram. We 
find no statistically significant (> la) correlation between 
the IR excess and M^s using the Spearman rank correla- 
tion test. This lack of correlation is also seen for the more 
homogeneous IC 348 core sample. (The exception is the 
M6.5— M7.5 bin, which exhibits a correlation of modest 
statistical significance. But this is driven by the one Tau- 
rus object with a very large E(Ks — L')q. When this ob- 
ject is excluded, the correlation disappears.) If we adopt 
an age for the substellar population based on that esti- 
mated for the stars (see § 2.1), our findings suggests that 
the inner regions (^0.1 AU) of young brown dwarf disks 
do not evolve substantially over the first ~3 Myr. This 
timescale is in accord with studies of disks around T Tauri 
stars (e.g. Strom et al. 1989; Skrutskie et al. 1990). 
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4. ANALYSIS OF DISK PROPERTIES 

4.1. Disk Models 

We use simple models of circumstellar disks to exam- 
ine the observed infrared excesses. Our goal here is not 
to deduce detailed physical properties, but to understand 
the general characteristics and trends of our observations. 
For simplicity, we assume that the disks are optically thick 
and geometrically flat. Similar models have been widely 
used to understand the emission from disks aroimd the 
higher mass T Tauri stars (e.g. Lynden-Bell & Pringle 
1974; Adams et al. 1987; Lada & Adams 1992; Hartmann 
1998, and references therein). The assumption of a flat 
disk is appropriate for the inner disk region from which 
the Ks — L' emission originates. Although in reality disks 
probably flare significantly at larger radii, the impact of 
the flaring on the SED is much more important at mid-IR 
(-^10 jjLm) and longer wavelengths (e.g. Kenyon & Hart- 
mann 1987; Wood et al. 2002). Indeed there are some in- 
dications from mid-IR measurcrncnts that flat disks might 
be more appropriate for some young brown dwarfs (Natta 
et al. 2002). 

For both active accretion and passive reprocessing of 
stellar irradiation, the disk temperature profile follows 
T oc r"'^/'', where r is the radial position in the disk. 
The normalization of the temperature profile is determined 
from the total luminosity that is produced by reprocessing 
and accretion. (The heating of the star by radiation from 
disk is ignored.) For a disk that extends from an infinite 
outer radius to an inner radius the luminosity from 
reprocessing is 



LfQ-p — 0.25 



(4) 



where L* and i?* are the luminosity and radius of the cen- 
tral source, respectively (Adams & Shu 1986). For the 
same disk, the luminosity from accretion is given by 



Lac 



(5) 



where M* is the mass of the central source and M is the 
disk accretion rate. Half of the total potential energy of 
the accreting material is radiated away, with the other half 
stored as rotational energy. In principle, the accretion lu- 
minosity can become arbitrarily large, since it scales with 
M, in contrast to the reprocessing luminosity which has 
a maximum of 0.25L*. In practice, for typical T Tauri 
stars, the luminosity in excess of the stellar photospheres 
is « 0.25L*, i.e., about the maximum expected from a flat 
reprocessing disk alone (Hartmann 1998). 

Given the contributions from reprocessing and accre- 
tion, we compute the resulting disk SED by summing over 
annuli, assuming each annulus radiates as a blackbody. To 
compute broad-band magnitudes from the resulting model 
SEDs, we use transmission proflles which account for the 
combined spectral response of the detector, filter, and at- 
mosphere. We obtained JHKs proflles from 2MASS, and 
L'M' ones for the MKO system from UKIRT (S. Leggett, 
priv. comm.). In computing the emergent flux from the 
star-|-disk spectrum, we account for the non-blackbody na- 
ture of the stellar spectra, a small, but noticeable, effect for 
cool M-type photospheres. We adjust the stellar spectra in 
the models to match the IR color locus for fleld M dwarfs 



described in § 2.3. In doing so, we make use of bolomet- 
ric corrections computed by Leggett et al. (2000, 2001). 
Although we have ignored, for simplicity, the true non- 
blackbody nature of the disk emission (e.g. Calvet et al. 
1992), we note that the results of these simple models agree 
well with those of Calvet et al. as presented in Meyer et al. 
(1997) for the low accretion rates considered here. 

Table 3 tabulates the maximum possible excess from our 
reprocessing disks, i.e., for a face-on orientation with no 
inner hole. Figure 6 shows the resulting JHKsL' color- 
color diagram for stars with reprocessing disks for a range 
of central stars, viewing angles, and hole sizes. The re- 
sulting models all have roughly the same color-color slope. 
The plotted stellar loci for giant and dwarf stars are from 
Besscll & Brett (1988), converted to the 2MASS system. 
We use 2nd order polynomial fits to represent the mean 
colors of Ml to M9.5 dwarfs, based on Leggett (1992), 
Leggett et al. (1998), and Leggett et al. (2002) data; the 
rms scatter about the fit is consistent with the measure- 
ment errors. We adopt the reddening law from Mathis 
(2000) for Ry = 3.1. Notice that (1) the disks around ear- 
lier type (hotter, more luminous) objects generate much 
larger IR excesses than those around later type (cooler, 
less luminous) objects; (2) mid-M type objects (near the 
stellar/substellar boundary) with disks are predicted to be 
bluer in J — H than either early-M types (representative 
of T Tauri stars) or late-M (substellar) objects with disks; 
and (3) models with significant holes generate color ex- 
cesses in only Ks — L', because the inner disk regions are 
relatively cool. 

4.2. Evidence for Inner Disk Holes 

L'-band observations are sensitive to emission from the 
hot, inner portions of circumstellar disks. For the range 
of Teff considered here, the emission comes from a region 
within a few radii of the central source (<0.1 AU). Hence, 
the presence of inner holes in the disks will have a substan- 
tial impact on the Ks — L' emission. Disks with inner holes 
of a few stellar radii have been inferred to exist around the 
higher mass T Tauri stars (e.g. Kenyon et al. 1996; Meyer 
et al. 1997) and are believed to originate from the inter- 
action of the stellar magnetic field with the accretion disk 
(e.g. Bertout et al. 1988; Konigl 1991; Shu et al. 1994). 
Modeling and interpreting the disk emission from a spe- 
cific object depends on the viewing angle to the observer, 
which is unknown. This limitation can be overcome with 
a large unbiased sample of objects, as is the case for our 
study — meaningful constraints on the inner holes can 
then be found from (1) the maximum observed IR excess 
and (2) the observed distribution of IR excesses. 

Figure 7 illustrates the effect of inner holes on passive 
disk models viewed face-on. Holes of >4i?* are sufficient 
to quench any Ks — L' excess, especially for later-type 
(cooler) objects. One interesting aspect of Figure 7 is that 
none of our targets are observed to have an IR excess as 
large as the model without an inner hole {Rin = R*)- 
The observed upper envelope is more consistent with the 
Rin ~ 2 — 3i?j, models, suggesting inner holes are com- 
mon. Note that for a specified hole size, the distribution 
of excesses should peak near the maximum excess based on 
simple geometrical considerations. The probability distri- 
bution of viewing angles follows P{9) = sin 6. Although 
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disks seen nearly face-on {9 ~ 0°) are rare, the decrease of 
disk flux with cos 6* compensates for this. The net result 
is that the color distribution in magnitudes for a random 
population of disks is strongly peaked towards the maxi- 
mum excesses (e.g. Kenyon & Hartmann 1990). 

To quantitatively examine the possibility of inner holes, 
we compare the observed E{Ks — L')o color distribution 
against Monte Carlo realizations of star+disk models 
viewed at randomly selected viewing angles. We first con- 
struct a scries of passive disk models with a fixed inner 
radius. The amplitude of the IR excess depends on the 
viewing angle and the T^g of the central source: hotter 
stars with disks can produce larger IR excesses than cooler 
stars with disks. However for the range of temperatures 
considered here, the Tcs dependence is removed when the 
model excesses are normalized against the maximum pos- 
sible excess (i.e., ^ = 0°). This normalization allows us to 
examine models which depend only on the viewing angle, 
which are needed to compare with the Monte Carlo simu- 
lations. Thus, for each object, we compute the normalized 
excess, defined as 



£{Ks-L')o 



Tpohs 



(6) 



max(i;f°l^|,) 

where the numerator is the observed {Ks — L')q excess and 
the denominator is the Kg — L' excess from a face-on pas- 
sive disk model computed for the Teff corresponding to 

the object's spectral type. £{Ks~L')q is simply the ra- 
tio of the observed excess to the maximum possible model 
excess. 

We exclude objects with small observed E{Ks — L')q 
from this analysis; these objects would correspond to 
model disks viewed nearly edgc-on. In such cases, the 
disk probably obscures the central source, and hence these 
objects are missing from the original imaging surveys. Ob- 
jects in our sample with small E{Ks — L')o are probably 
not those with disks viewed nearly edge-on, but rather ob- 
jects without significant disks. We choose a cutoff value 
of ^ = 80°, corresponding to a 100 AU fiared disk with an 
outer height of 17 AU; this agrees well with results from 
HST imaging of the edge-on Taurus object HH 30 (Bur- 
rows et al. 1996). 

For a given inner hole radius, we compare the observed 
£{Ks — L')q distribution with that expected from a Monte 
Carlo set of disks chosen with random viewing angles. The 
results are shown in Figure 8 for models with inner disk 
radii of 1, 2, and 3i?^. Models without inner holes predict 
many more objects with large values of £{Ks — L')o than 
observed; this is equivalent to the findings in Figure 7, 
where there are no objects with excesses as large as the 
faco-on models without holes. On the other hand, models 
with inner radii of 3i?* predict much smaller excesses than 
observed. Only models with inner radii of w 2Ri, agree 
well with the data. This is confirmed using the K-S test: 
for a hole radius of 2.2Ri,, the K-S test indicates a 82% 
chance that the observed and the predicted £{Ks — L')o 
distributions originate from the same population. This is 
much better concordance than for the li?* and 3i?* hole 
models; the K-S test indicates a > 1000 x lower proba- 
bility of these models being drawn from the same parent 
population as the observations. 

To summarize, we find the observed maximum ampli- 
tude and color distribution of the {Ks — L')o excesses point 



to disks with characteristic inner radii of « 2i?*. Note that 
our analysis based on passive disks provides a lower esti- 
mate of the hole size: larger holes could be accommodated 
if the luminosity from accretion is substantial (see § 4.3) 
or if the vertical flaring of the disks is unexpectedly large. 
Mid-IR observations of young brown dwarfs in p Oph may 
also suggest the presence of inner holes (Natta et al. 2002). 
Finally, as illustrated in Figure 7, holes of > 2Ri, could ac- 
count for the small/negligible IR excesses observed for the 
M9 sources in our sample. Longer wavelength data are 
need to conclusively determine whether these very cool 
objects have disks. 

4.3. Constraints on Accretion Rates 

The widespread presence of Ha emission among the 
sample (§ 3.2) indicates that accretion is ongoing for most 
objects. To roughly constrain the accretion rate, we rear- 
range equations (4) and (5) to compare the relative con- 
tributions of reprocessing and accretion luminosity for the 
case of a simple flat accretion disk: 



-'rep 




, , , (7) 

^0.02^0^ \0.6RqJ ^ ' 

where we have used values for the mass, luminosity and ra- 
dius representative of our sample. As discussed in the pre- 
vious section, the magnitude of the observed IR excesses 
can be generated by a reprocessing disk alone. There is no 
need to invoke significant luminosity from mass accretion, 
which suggests accretion rates are < 10~^ Mq yr~^. Such 
rates arc comparable to the lowest accretion rates mea- 
sured for classical T Tauri stars (CTTS); typical CTTS 
rates are around 10"'' to 10~* Mq yr~^ (e.g. Valenti et al. 
1993; GuUbring et al. 1998). Such low rates also indicate 
that young brown dwarfs accumulate at most ~10% of 
their final mass via disk accretion, if the current rates are 
representative of their time-averaged accretion history. 

For one of our targets, the M6 object V410 Anon 13 in 
Taurus, MuzeroUe et al. (2000) find its Ha line profile is 
well-fitted using magnetospheric accretion models in which 
the magnetospherc extends to 2.2— 3ii*, in accord with our 
analysis of the {Kg — L')o excesses in the previous section. 
They derive an accretion rate of ~ 5 x 10~^^ Mq yr~^, 
well below our estimated upper limit. More such direct 
measurements for young brown dwarfs are sorely needed 
to determine the typical accretion rate and its dispersion. 

As mentioned before, there is a degeneracy in determin- 
ing the inner hole size and the accretion rate based on the 
IR data alone. Pure reprocessing disks with Ki2Ri, inner 
holes can account for the observed IR excesses (§ 4.2). 
But disks with larger inner holes can also agree with 
the data if the accretion rates are significant. For in- 
stance, disk models with 3iZ^, holes and accretion rates of 
^10^" Mq yr^^ can generate comparable excesses. How- 
ever, such accretion rates would disagree with Ha measure- 
ments of young brown dwarfs by Muzerolle et al. (2000, 
2002), which include a few objects common to our V- 
band sample. Even larger inner holes of 4— Si?^, would 
make reprocessing nearly irrelevant (e.g.. Figure 7), and 
the accretion rates required to generate the observed ex- 
cesses (>10~^ Mq yr~^) would be even more discrepant. 
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Thus it is unlikely that the characteristic sizes of the inner 
holes are much larger than Ri2i?^. 

5. POTENTIAL SYSTEMATIC EFFECTS 

Our identification and analysis of IR excesses depend on 
several measurements which may be prone to systematic 
errors. Here, we find that they should not have a signifi- 
cant impact. 

5.1. Extinctions and Spectral Types 

Measuring the IR excesses depends on the adopted spec- 
tral types and extinctions. As a check on these, we ex- 
amine the subset of our IC 348 core sample which has 
been studied by more than one group. Twelve of the ob- 
jects, about half the sample, have been classified using 
both IR narrow-band spectrophotometry (Najita ct al. 
2000) and optical/far-red spectroscopy (Luhman 1999); 
these twelve are roughly distributed across our sample's 
full range of spectral types. The ii'-band extinctions mea- 
sured by the two studies agree well, with a median differ- 
ence of 0.03 mag and an rms of 0.12 mag. This concurs 
with the analysis of Najita et al. for their entire IC 348 
sample. Note that since we are measuring IR color ex- 
cesses, we are somewhat less sensitive to errors in Ak than 
if we were simply measuring bandpass magnitudes. 

As for the spectral types, the Luhman spectral types 
tend to be «1 subclass earlier than those of Najita et al. 
This is expected, given the spectral type calibration chosen 
by Najita et al. (see their Figure 7, bottom panel). If the 
Luhman (1999) types are adopted, the net result is that 
the IR excesses we measure would be larger, strengthening 
our finding that IR excesses are common in this sample. 
Of course, the mass estimates would be increased, pushing 
some objects from brown dwarf to VLM stars, but this 
would not change our findings about the high frequency 
of IR excesses in this very low mass sample. Likewise, 
our analysis of the inner hole sizes (§ 4.2) would not be 
impacted much, since the maximum E{Ks — L')o for the 
models has only a weak dependence on spectral type for 
M5 to M9 objects (Figure 7). Hence, the uncertainties in 
the spectral types and extinctions do not seriously affect 
our analysis. 

5.2. Contamination from Disk Emission 

Disk emission can in principal lead to systematic errors 
in the derived spectral types and extinctions. This is un- 
likely to be an issue for the Taurus and IC 348 outer sam- 
ples, since the spectral typing was done at optical wave- 
lengths (7-band) where any disk emission should be negli- 
gible. However, disk emission contamination is potentially 
more important for the IC 348 core sample, whose spec- 
tral types are based on NICMOS 1.6—2.2 /im narrow-band 
photometry measuring the 1.9 /^m H2O absorption band. 
In a standard passive and/or accretion disk where the flux 
excess follows the form AF\ cx X~^/^, disk emission would 
lead to weaker H2O absorption, and hence earlier spectral 
types. 

We use our passive disk models to quantify the impact 
of disk emission on the NICMOS spectral types. We ex- 
plore the worst-casc! scenario, in which the disks do not 
have inner holes and hence the disk emission is maximized. 
Our approach is as follows: for each object, we construct 



passive disk models based on the measured spectral type 
and spanning a wide range of viewing angle 6. We use 
the observed E[Ks — L')q to assign a corresponding value 
of 9 and then determine the amount of disk excess in the 
NICMOS narrow-band filters. Employing the Najita et al. 
calibration of NICMOS photometry to spectral type, we 
then determine the change in typing. (For a fully con- 
sistent analysis, we would do this iteratively, using the 
inferred error in spectral type to recompute the observed 
E{Ks — L')a, and then repeating the above process to re- 
determine the systematic error in spectral type. In prac- 
tice, this is unnecessary since the disk contamination is 
small.) 

For the IC 348 core sample, we find the median effect of 
contamination by disk emission is 0.2 subclasses, i.e., the 
objects should be classified as slightly later types. The me- 
dian excess in the ii'-band photometry is 0.07 mags. As 
expected, the objects with the most significant disk con- 
tamination are those with the earliest spectral types and 
largest E{Ks — L')q: here, the effect of disk emission can 
be as large as w 0.5 — 0.7 subclasses, producing a i^-band 
excess of ~ 0.2 — 0.3 mags. Even in these cases, the conse- 
quences arc minor: the effect of a later spectral type is to 
reduce the measured E{Ks — L')q. However, the actual im- 
pact is small since the photospheric Ks — L' color changes 
slowly with spectral type for late M-dwarfs. Therefore, we 
conclude that the contamination from disk emission is not 
significant. In a similar vein, the impact of disk emission 
on the extinctions derived by Najita et al. (2000) is small. 

Again, since these calculations use models without in- 
ner holes, the above results present the worst case scenario 
for disk contamination. In the event that the disks have 
inner holes of radii > 2i?*, as we infer in § 4.2, the effect of 
disk emission on the NICMOS-derived spectral types and 
extinctions is negligible. 

5.3. Surface Gravity 

We have used the colors of field M dwarfs in order 
to identify which of our targets have intrinsic IR ex- 
cesses. However, our targets are pre-main sequence ob- 
jects, and hence will have lower surface gravities than the 
field M dwarfs. The sense of this effect is likely that at ob- 
jects with lower gravities have bluer Ks — L' colors. This 
is seen in Figure 6, where the Ks — L' colors of giant stars 
(spectral types GO HI to M5 HI) are comparable to or bluer 
than dwarf stars (AO V to M9.5 V). Pre-main sequence 
stars are intermediate in surface gravity, and hence we ex- 
pect their Ks — L' colors to be no redder than the dwarf 
stars. If they are bluer than dwarfs, then we would infer 
that the IR excesses are larger and even more common. 

Given the current uncertainties in the theoretical at- 
mospheres, we chose not to use the models to measure 
E{Ks-~ L')q in § 3, but instead used the empirical field 
dwarf locus. Here, we examine the Baraffe et al. (1998, 
2002) models to explore the relative shift in color due to 
gravity effects. Figure 9 plots models with ages of 1 Myr 
and 1 Gyr: there is a difference of a factor of w30 in 
the surface gravities at fixed mass. The T^g for the older 
(1 Gyr) models have been converted to spectral type using 
results for field M dwarfs from Leggett et al. (2000, 2001) 
and Tsuji et al. (1996), fitted to a hnear relation: 

Teff{K) = 3851 - 185.3 x SpT, (8) 
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where SpT is the M spectral subclass. For the younger 
(1 Myr) models we use either this scale or the aforemen- 
tioned Luhman (1999) scale (equation (1)). Both the 
"clear" and the "dusty" atmospheric models predict that 
younger, lower gravity objects have bluer colors at fixed 
spectral type. The amplitude of the difference varies with 
the choice of model and temperature scale, but generally 
lower gravity objects are predicted to be 0.1 mag bluer in 
Ks — L' at fixed spectral type. 

Perhaps the most convincing evidence that gravity ef- 
fects do not significantly affect the photospheric colors can 
be seen in the data themselves. The loci of the observed 
and dereddened Kg — L' colors in Figure 2 both show a 
lower envelope which agrees well with the empirical dwarf 
locus, to better than < 0.05 mag. This level of difiierence 
has little bearing on our analysis.^ 

5.4. Sample Selection Biases 

Biases in the sample selection could impact our con- 
clusions. For example, if brown dwarfs without disks are 
systematically excluded from our sample, then our esti- 
mate of the disk-bearing fraction will be incorrect. The 
effect of extinction on the sample selection is probably the 
primary concern. Objects with edge-on disks will be heav- 
ily extincted and are almost certainly missing from the 
optical/IR imaging surveys used to construct our sample. 
(Indeed, this bias is also present in most previous T Tauri 
disk studies using broad-band photometry.) Viewing an- 
gles of greater than 80°, which is the cutoff value we adopt 
in § 4.2 for obscured edge-on disks, account for 17% of a 
random distribution of disk inclinations, so the disk frac- 
tion for our sample would be underestimated by only a 
smaU amount (w3%). Also, including such objects would 
only act to increase the measured disk fraction, which we 
already find is very high. 

One common approach to minimizing selection biases is 
to use a complete sample, where all the objects in a given 
area on the sky are identified to a well-defined limit. The 
Najita et al. (2000) survey is the source of our IC 348 core 
sample; this survey is if-band selected and has a well- 
defined completeness limit for spectral classification. We 
have i'-band observations for nearly every cluster member 
from their survey in the M5.7 to M9.4 range; only 2 objects 
out of 26 were not observed. Also, most of our targets are 
several magnitudes brighter in if-band than the survey 
limit. Therefore, we expect that our IC 348 core sample 
is essentially complete, and that selection biases due to 
extinction are negligible. The trends we identified in § 3 
for our entire L'-band sample are also found if we consider 
only the IC 348 core sample. This suggests that selection 
biases are probably not a significant effect when analyzing 
the L'-band sample as a whole. 

The three coolest/lowest mass objects, all with spec- 
tral type «M9, have negligible IR excesses. This might 
be due to a selection effect if objects with larger excesses 
tended to be preferentially extincted and hence missing 

* One interesting feature is that the dusty models predict fairly constant K — L' colors for objects with mid to late-M type objects, espe- 
cially when using the Luhman (1999) temperature scale. This is inconsistent with the observations. For field M6— M9 dwarfs, the models are 
Ri 0.1 — 0.2 mag bluer than the observed colors (e.g., the solid line plotted in Figure 2). Likewise, the dusty models have bluer Ks — L' colors 
than the observed lower envelope for our sample of young M dwarfs (Figure 2). However, our sample contains only a few young M9 objects, 
so the observations would allow for constant K — L' colors if this effect in fact begins around M8, instead of at M6 as currently predicted by 
the dusty models. If this is the case, then the M9 objects would have actually small IR excesses, whereas they are currently inferred to have 
almost none. 



from the sample. But there is no evidence for this ef- 
fect: for the M5.7 to M7.7 objects, no statistically signif- 
icant (> la) correlation exists between E{Ks — L')o and 
Ak- Also, even the M9 objects are 1.5 mags above the 
K-hand survey limit of Najita et al. Thus, it is unlikely 
there are M9 objects with large excesses which are missing 
from our sample because of large extinctions. However, a 
much larger sample of M9 (and cooler) sources are needed 
to accrurately assess the disk fraction of such low mass 
objects (< 15-20Mj„p). 

5.5. Pre-Main Sequence Temperature Scale and Mass 

Estimates 

When analyzing the trends in the disk properties (§ 3), 
we used the observable quantities of spectral type and 
Mks as surrogates for the relative masses and ages of the 
targets. Preferably, one would like to use absolute values 
for masses and ages, as obtained from theoretical models. 
However, an examination of Figure 1 shows the problem: 
different models give very different results. Also, the cur- 
rent models do not always cover the observed location of 
the young brown dwarfs in the plot. (This is true using ei- 
ther K-haiiA absolute magnitude or bolometric magnitude 
as the ordinate.) 

An important caveat is the choice of the conversion be- 
tween Toff and spectral type. We have used the Luhman 
(1999) scale throughout. If we had used a scale appro- 
priate for field M dwarfs (equation (8)) then evolutionary 
models of a given temperature would correspond to earlier 
spectral types — the net result would be much lower mass 
estimates for our targets, with masses of <10 Mjup for the 
coolest ones. Hence our choice of the Luhman (1999) scale 
is a conservative one. Similarly, the disk modeling results 
in Table 3 would change by < 0.05 mag with a dwarf tem- 
perature scale, and our other analyses using disk models 
in § 4 would also be unaffected. 

We reiterate that our targets are selected based on their 
spectral types. The resulting sample is significantly cooler 
than previous L-band surveys, which were sensitive down 
to spectral types of ^M5 (Kenyon & Hartmann 1995; 
Haisch et al. 2001a), and hence comprise the lowest mass 
sample of young objects studied to date at these wave- 
lengths, independent of the choice of temperature scale. 

6. DISCUSSION 

6.1. Disk Frequency of Young Stars and Brown Dwarfs 

In our sample of brown dwarfs and VLM stars, we 
find 31 out of 38 objects possess intrinsic IR excesses 
(82% ± 15%, where we quote only Poisson statistical er- 
rors in this section in order to compare with previously 
published results). This disk fraction is comparable to, or 
exceeds, the disk fraction for the higher mass stars in the 
same star-forming regions. Using multi-band photometry 
and spectral types, Strom et al. (1989) found 47 out of 
83 (57% ± 8%i) T Tauri stars in Taurus have significant 
if- band excesses {AK > 0.1 mag). (See also Kenyon & 
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Hartmann 1995.) Haisch et al. (2001a) found a similar 
disk fraction of 65% ± 8% for stars in IC 348. Some cau- 
tion is in order when comparing our IC 348 results with 
those of Haisch et al. The latter used only multi-band pho- 
tometry without spectral type information, and hence, as 
we discuss below, their disk fraction may be somewhat 
underestimated. 

The most common method for measuring the disk frac- 
tion of young stars uses color-color diagrams based on 
JHK, or preferably JHKL, colors. Objects with disks 
are identified as those having IR colors distinct from red- 
dened dwarf and giant stars. This method has been used 
extensively to study the frequency, properties, and evo- 
lution of circumstellar disks around T Tauri stars (e.g. 
Lada & Adams 1992: Kcnyon & Hartmann 1995; Lada & 
Lada 1995; Meyer et al. 1997; Haisch et al. 2001b). The 
method is appealing since only photometry is used, with- 
out need for more time-consuming spectroscopy. More- 
over, the method has been demonstrated to be effective for 
T Tauri stars, e.g., Hillenbrand et al. (1998) and Haisch 
et al. (2001c) show that most TTS in Taurus with multi- 
wavelength evidence for disks can be identified from JHK 
colors alone, and all can be found with JHKL colors. Here 
we consider the value of this method for studying disks 
around young brown dwarfs. 

Figure 10 illustrates the commonly used IR color-color 
analysis applied to our sample. Prom this diagram, one 
would identify only 11 out of 36 objects as having IR ex- 
cesses. (Two objects do not have 2MASS J-band pho- 
tometry.) However, our analysis which incorporates the 
objects' spectral types and extinctions shows that in fact 
many more objects (31 out of 36) have IR excesses. The 
majority of sources with disks arc missed in the color-color 
diagram because their IR excesses are modest. In fact, an 
analysis using only JHKg colors, without i'-band data, 
would find only 2 out of 36 objects with IR excesses, i.e., 
essentially missing all the objects with disks. These ob- 
servational results confirm model predictions of very small 
JHK excesses from brown dwarf disks (e.g., Natta & Testi 
2001, this work). 

There are two physical reasons why using only IR col- 
ors works poorly. Both lead to decreased contrast between 
brown dwarfs and their disks. (1) Because brown dwarfs 
are less luminous than the higher mass T Tauri stars, disks 
around brown dwarfs will be cooler and less luminous; 
therefore, the corresponding IR excesses will be smaller 
(e.g., Table 3). (2) T Tauri stars span a hmited range 
in photospheric Ks — L' color, and hence even a modest 
disk excess will produce IR colors readily distinguishable 
in color-color diagrams. In contrast, young brown dwarfs, 
with spectral types of M6 and later, have redder photo- 
spheres and span a larger range in intrinsic Ks — L' color 
(e.g., Figure 6). Hence, the modest IR excesses produced 
by their disks will be harder to distinguish than for the 
case of T Tauri stars. 

This latter item raises the fact that color-color analyses 
inevitably underestimate the disk fraction. Near-IR imag- 
ing samples will include objects with a range of spectral 
types. Earlier-type (hotter) stars with small IR excesses 
will not be as readily identified as disk-bearing in color- 
color diagrams compared to stars with larger excesses or 
later-type (cooler) objects which lie close to the reddening 



boundary. This effect leads to an underestimate of the disk 
fraction. It is unlikely to be an significant error in surveys 
of T Tauri or Herbig Ae/Be stars, since they have a small 
range in intrinsic IR colors. It will be a more serious effect 
for young brown dwarfs. 

In a similar vein, the disk fractions inferred from color- 
color analyses are very sensitive to the assumed boundary 
of the reddened stellar locus (the rightmost line in Fig- 
ure 10). The boundary is determined by the colors of the 
latest spectral type objects. However, since spectroscopy 
is generally not available in these studies, this boundary 
must be estimated, which can produce either an under- 
estimate or overestimate of the true disk fraction. Since 
nearly all such surveys are magnitude-limited, this sys- 
tematic error can be aggravated by the variable amounts 
of extinction to different sample members and hence an 
ill-defined spectral type boundary. 

Our approach of examining a sample of objects with 
spectral type determinations is far more sensitive to disk 
emission than searches based on JHK{L) photometry 
alone. Indeed, our results suggest that both L'-band data 
and spectral types are required in order for an accurate 
ground-based census of disks around young brown dwarfs. 

6.2. Comparison with Near-IR Surveys of Trapezium and 

a Ori 

Two recent surveys have addressed the topic of disks 
around young VLM stars and brown dwarfs using large 
samples in Orion. Here we compare their findings with 
our own results. 

Muench et al. (2001) have studied the disk fraction 
for young brown dwarfs in the Trapezium cluster using 
JHK imaging. They select a large sample (109 objects) of 
brown dwarf candidates identified from their {J — H,H) 
color-magnitude diagram in conjunction with a 1 Myr 
theoretical model from Baraffe et al. (1998). Then they 
use JHK color-color diagrams to infer a disk fraction of 
~ 65% ± 15%. Many of their objects with IR excesses 
clearly possess disks, because they are coincident with 
"proplyd" sources seen in HST optical images. 

Interpreting their results is complicated by the uncer- 
tain masses of their candidates. Contamination from back- 
ground stars is predicted to be low for the Muench et al. 
magnitude limits (e.g. Hillenbrand & Carpenter 2000), so 
most of the objects probably are cluster members. But 
brown dwarf selection from the IR color-magnitude dia- 
gram is clearly not robust: about 15% of their candidates 
have JHK colors consistent with reddened stars of spectral 
types MO or earlier. Moreover, there is a puzzling inconsis- 
tency in the photometric selection when compared to the 
available spectroscopy. Ten sources in the Muench et al. 
survey have spectral types of M6 or later as measured by 
Hillenbrand (1997), indicating that they are near or below 
the substellar limit. However, all of these are brighter than 
the 1 Myr model predictions for a 0.08 Mq object, some 
by several magnitudes. The evolutionary models would 
thus indicate that these objects have stellar masses based 
on their brightnesses and are not brown dwarfs. Finally, 
the Muench et al. selection approach essentially uses un- 
extincted J-band magnitudes, and therefore depends on 
the assumed age of the Trapezium cluster. Objects older 
than 1 Myr, and hence higher mass, might be included in 
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the sample. For example, for an assumed age of 3 Myr, 
the brightest candidates would correspond to a mass of 

0. 15 Mq based on the Baraffe et al. (1998) models. Spec- 
troscopic follow-up will be important for refining the mass 
estimates of their brown dwarf candidates. 

The location of the Muench et al. Trapezium brown 
dwarf candidates in the JHK color-color diagram is in- 
triguing in light of our modeling results. (1) Many of the 
objects exhibit H—K excesses of several tenths of a mag, 
substantially larger than the maximum excess predicted 
by simple flat disk models (Table 3). One possible expla- 
nation is that these are flared disks seen close to edge-on; 
models which include the effects of scattering show that 
such disks tend to have redder colors than simple flat re- 
processing disks (Wood et al. 2002). If this is the case, the 
central sources will be more heavily extincted, and since 
the target selection was based on magnitudes, will tend to 
have higher masses. Another possibility is that accretion 
luminosity is a significant component of the disk emission, 

1. e., M > 10"^ Mq yr-i (equation (7)). In this case, the 
implication would be much higher accretion rates for the 
Trapezium sample compared to the Taurus and IC 348 
objects in our sample, perhaps a reflection of Trapezium's 
younger age and/or higher stellar density. (2) Most of the 
Trapezium IR excess sources lie below the CTTS locus of 
Meyer et al. (1997), namely they have bluer (J—H) colors. 
Simple disk models can explain some of the objects as be- 
ing mid-M type objects with disks, which will have bluer 
J—H colors than the early-M type objects used to define 
the CTTS locus (e.g.. Figure 6). However, these models 
cannot explain sources with very blue colors {J—H < 0.6 
mag). Muench et al. point out that half of these are asso- 
ciated with proplyds and hence are gc;nuine disk-bearing 
sources; perhaps their blue colors arise from scattering. 
The nature of the remaining blue sources remains unre- 
solved. 

Oliveira et al. (2002) have studied the disk fraction for 
low mass stars and a few brown dwarf candidates in the 
£7 Ori cluster. For a sample of 34 cluster members se- 
lected on the basis of optical photometry and Li absorp- 
tion, they found at most 2 out of 34 (6% ± 4%) objects 
had IR excesses, based on JHK colors. Adopting a clus- 
ter age of ~2-7 Myr, Oliveira et al. suggest that disks 
around very low mass objects can dissipate within a few 
Myr. These conclusions are cast into doubt by our results. 
We have shown that JHK colors arc; poor diagnostics of 
disks around VLM stars and brown dwarfs. Longer wave- 
length observations, at least in the L-band, are needed 
to measure the true disk fraction and its implications for 
the evolutionary time scale of disks around very low mass 
objects. 

6.3. Properties of Young Brown Dwarfs and their 

Circum(suh) stellar Disks 

Our results suggest that much of the observational 
paradigm developed for disks around T Tauri stars can 
be extended to disks around young brown dwarfs. Most 
of our sample show both IR excesses and strong Ha emis- 
sion, like the higher mass classical T Tauri stars. Some 

^ This evidence for strong magnetic activity among young brown dwarfs contrasts strongly with the evidence for decreased magnetic activity 
among dwarfs at the bottom of the main sequence (Gizis et al. 2000). This is likely due to the higher temperatures of young brown dwarfs (i.e., 
spectral types of mid-M and later), which are high enough to support a partially ionized atmosphere, and therefore, surface magnetic fields, 
chromospheres, and coronae (e.g. Preibisch & Zinnecker 2002). 



objects have little/no IR excesses and little Ha emission, 
analogous to weak-line T Tauri stars. A few have IR ex- 
cesses but only weak Ha emission; such a phenomenon is 
also seen among weak-line T Tauri stars (e.g. Strom et al. 
1989; Skrutskie et al. 1990; Preibisch & Zinnecker 2002). 
In addition, the widespread presence of Ha emission sug- 
gests that most brown dwarfs are accreting, although the 
inferred rates are much lower than classical T Tauri stars. 

We have also found evidence that young brown dwarfs, 
like their more massive counterparts, the T Tauri stars, 
possess circumstellar disks with inner holes (§ 4.2). What 
is the origin of these holes? In the case of T Tauri stars, in- 
ner holes are believed to be created through the truncation 
of the disk by strong, closed stellar magnetic fields; disk 
matter may reach the star by accretion along these stellar 
magnetic field lines (e.g. Bertout et al. 1988; Konigl 1991; 
Hartmann et al. 1994; Shu et al. 1994). A similar situation 
may be relevant to young brown dwarfs. Like their more 
massive counterparts, young brown dwarfs also show signs 
of surface magnetic activity. For example, young brown 
dwarfs have detectable X-ray emission that is character- 
ized by plasma temperatures of ~ 1—2 keV and X-ray lumi- 
nosities that are ^ 10-"* to 10''^ times the bolometric lu- 
minosity. The X-ray emission also displays moderate time 
variability and rapid flaring (Imanishi et al. 2001; Preibisch 
& Zinnecker 2001, 2002). These similarities strongly sug- 
gest that young brown dwarfs have significant magnetic 
fields.5 

The likely disk accretion rates of the brown dwarfs in 
our survey are expected to be small (§ 4.3) and only mod- 
est magnetic fields are needed to produce the inner holes 
with radii of 2—3 i?*. For example, with typical param- 
eters of M, = 0.06 Mq, i?, = 0.6 i?0, and a disk accre- 
tion rate of 10~^° Mq yr~^, the required field strengths 
are ~ 110 - 220 G if the X-wind theory of Shu et al. 
(1994, equation 2.6b with ax ~ 1) is applied directly to 
this case (sec also Ostriker & Shu 1995 who recommend 
ax = 0.923). This estimate assumes that the magneto- 
sphere is only partially filled, with the accretion rate onto 
the star restricted to stellar loops that correspond to ~10% 
of a filled magnetosphere (Calvet & GuUbring 1998; Muze- 
roUe et al. 2000). For comparison, the formulation given by 
Konigl (1991) when applied directly to this case requires 
field strengths of ~ 450 - 900 G. These field strengths 
are quite modest, less than those measured for pre-main- 
sequencc stars (Gucnther et al. 1999; Johns-KruU et al. 
1999a,b) and dMe stars (e.g. Johns-KruU & Valenti 1996). 

It is interesting that the temperature at an inner disk 
radius of 2i?* in our disk models is approximately 1200 K, 
very similar to the estimated temperature for the inner 
disk radii of T Tauri stars (e.g. Shu et al. 1994). At the 
densities of inner disks and temperatures of >;1000 K, ther- 
mal ionization (of Na and K) can maintain an ionization 
fraction that enables magnetic coupling between the star 
(or brown dwarf) and the disk (Umebayashi & Nakano 
1988; Fromang et al. 2002). Therefore, the inner regions 
of brown dwarf disks will be sufficiently ionized for mag- 
netic fields to couple to and truncate the disk at w 2J?j,. 
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6.4. Implications for Brown Dwarf Formation 

For a star like the Sun, significant disk accretion is an 
imperative because the angular momentum of the parent 
cloud core is typically too large for the star to accrete more 
than a few percent of its final mass through direct infall 
(e.g. Shu et al. 1993). If Sun-like stars and brown dwarfs 
form under similar initial conditions (e.g., they form from 
similar cloud cores), then the need for significant disk ac- 
cretion is much reduced since brown dwarfs may be able to 
accrete a substantial fraction of their final mass through 
direct infall. Our results indicate that if brown dwarfs form 
through such means then disks are produced. In such a 
scenario, an important question is what determines the fi- 
nal mass of the object, i.e., why a given core would create 
a brown dwarf rather than a Sun-like star. 

One physical effect which may address this concern is 
the premature truncation of accretion due to strong dy- 
namical interactions at early times. Reipurth & Clarke 
(2001) have proposed that brown dwarfs originate as stel- 
lar embryos in small multiple systems which are prema- 
turely ejected. These embryos accrete from a common in- 
falling envelope, but dynamical interactions among them 
lead to the preferential ejection of the lowest mass mem- 
bers. Here, the final masses of the ejected objects are 
established by the termination of accretion at the time of 
ejection. (Note that although their work focuses on brown 
dwarfs, there is no clear distinction across the substellar 
boundary and hence this scenario might also apply to VLM 
stars.) The process is purported to occur at very young 
ages (~10'' yr), due to the short crossing times of the mul- 
tiple systems. Preliminary studies of young star-forming 
regions do not find any obvious kinematic difference be- 
tween the brown dwarfs and T Tauri stars (Joergens & 
Guenther 2001), as might be expected from preferential 
ejection of the young brown dwarfs. However, Reipurth 
et al. (2001) suggest that the dynamical signature of this 
process may be subtle and largely independent of mass. 

An alternative class of formation scenarios invokes cir- 
cumstellar disks as the birthplace of substellar objects, 
perhaps through collisions of star-|-disk systems or intrin- 
sic instabilities in single star-|-disk systems (e.g. Lin et al. 
1998; Watkins et al. 1998a,b; Boss 1998). Specific predic- 
tions are lacking; e.g., for simulations of star-|-disk colli- 
sions, the masses of the resulting objects depend on the 
assumed star-l-disk properties and especially the orbital 
parameters. Nevertheless, dynamical interactions, either 
during the formation episode or afterwards, are needed 
to liberate the brown dwarfs into free-floating objects. 
In this respect, similarities exist with the aforementioned 
Reipurth & Clarke (2001) scenario. In a broader sense, all 
these scenarios have randomness as a central determinis- 
tic element, as opposed to, e.g., the initial conditions in 
pre-(sub)stellar molecular cores. 

Intuition suggests that ejection will be inherently hostile 
to disks around young substellar objects. Studies of the 
encounter of a star-l-disk system with another passing sys- 
tem find that such events are highly disruptive to the disks, 
typically truncating the disk sizes to half the periastron 
distance (Clarke & Pringle 1993; Hall et al. 1996). A sim- 
ilar outcome is seen by Armitage & Clarke (1997) for the 
ejection of a star-hdisk system from a small stellar cluster. 
These effects are likely to be more severe for disks around 



brown dwarfs given the shallower gravitational potential. 
Furthermore, such encounters are likely to trigger rapid ac- 
cretion (Bonnell & Bastien 1992). Hence, Reipurth (2000) 
and Reipurth & Clarke (2001) predict that young brown 
dwarfs are unlikely to have signatures of disks, or else that 
such signatures will have shorter lifetimes than for young 
stars. This is inconsistent with our results: we find that 
disks around young brown dwarfs are very common. More- 
over, the disks are contemporaneous with disks around the 
young stars in the same star- forming regions. Assuming 
the brown dwarfs and stars are roughly coeval, the life- 
times of brown dwarf disks are at least as long as those of 
stellar disks. 

One possible explanation for this discrepancy could be 
that most of the (diskless) ejcx",tc;d objects have dispersed 
to large distances from the young stars, which would oc- 
cur for ejection velocities of a few km s~^. Such objects 
would then be missing from existing optical/IR imaging 
surveys. The remaining brown dwarfs would be the ones 
which experienced only very gentle ejections, and hence 
their disks would be preserved. It is not clear if such a 
direct correspondence between ejection velocity and disk 
survival can be maintained. Moreover, in such a scenario 
the mass function measured for IC 348 would be strongly 
deficient in brown dwarfs compared to the field mass func- 
tion. This is not supported by the observed IMF (Luhman 
1999; Najita et al. 2000). 

Another possible explanation is that the post-ejection 
disks are truncated but they accrete at such low rates 
that they still retain enough disk material to produce an 
observable disk signature at ages of a few Myr. Detailed 
predictions for the effect of ejection on brown dwarf disks 
are not yet available, but simple estimates suggest that the 
surviving disks may be quite tenuous. We consider the case 
where the velocities of the ejected brown dwarfs are not 
very different than the stars, consistent with simulations of 
decaying few-body systems (Sterzik & Durisen 1998). For 
a 50 Mjup object, a typical velocity of ~2 km s~^ (.Jones 
& Herbig 1979; Hartmann et al. 1986) corresponds to the 
escape velocity at a distance of 10 AU. Adopting this es- 
timate as the truncation radius, we can ask whether the 
disks will have enough material to persist for a few Myr. 
The minimum total disk mass for the presence of optically 
thick L'-band emission is only ~10~^ Mq, which would 
quickly be accreted for canonical accretion rates of 10~^° 
to 10~" Mq yr^i. Instead assume that the initial disk was 
100 AU in size with a mass of 1% of the central object, 
similar to the case for T Tauri stars (Osterloh & Beck- 
with 1995). Assuming a surface density profile of r~'^/^ 
(Hayashi 1981), the truncated disk wiU have (10/100)^/2 
of the initial mass, or a total mass of about lO""* Mq 
for a 50 Mjup object. For the aforementioned accretion 
rates, the lifetime of these truncated disks will be about 
1-10 Myr, comparable to the ages of our targets. Hence, 
the disk mass budget is roughly compatible with the ejec- 
tion scenario. Note that this estimate is based on two 
assumptions: (1) the current accretion rates are roughly 
constant over the entire lifetime, and (2) the disk surface 
density profile is not perturbed significantly by the ejec- 
tion-htruncation event. Neither of these is likely to be 
realistic. 

Alternatively, our finding that brown dwarf disks are 
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as common as disks aroimd T Tauri stars is well- 
accommodated in perhaps the simplest picture of brown 
dwarf formation, namely the collapse of isolated, very low 
mass cores — a direct analog to the conventional picture 
of isolated low-mass star formation. In this case, the low 
disk accretion rates found for young brown dwarfs may 
have played a more central role in their formation. That 
is, brown dwarfs that coexist with stars in young clusters 
like IC 348 may have formed like their more massive coun- 
terparts but with much lower accretion rates in order for 
all the objects to have accumulated their final masses on 
roughly the same timescale. The idea that stars of dif- 
ferent masses result from different accretion rates is not 
new. For example, something similar is needed in order to 
explain the coexistence of high mass stars (e.g., >20 Mq) 
and low mass stars (e.g., < IMq) in young clusters (e.g. 
Shu et al. 1987; Myers & Fuller 1993). Since the main se- 
quence lifetime of high mass stars is only a few Myr, they 
must accrete their final mass on a timescale much shorter 
than this in order to coexist with much lower mass stars. 
The mass accretion rates required (~ 10"'' M© yr~^) are 
much larger than those typically ascribed to the young 
progenitors of solar-mass stars. Our results suggest that 
a similar scenario of mass-dependent accretion rates may 
apply down to substcllar masses. 

To summarize, the widespread presence of disks around 
young brown dwarfs at ages of a few Myr provides prima 
facie evidence of a similar origin for low-mass stars and 
brown dwarfs. The high disk frequency is more difficult 
to reconcile with scenarios involving disk-disk collisions 
and/or premature ejection. Important future observa- 
tional goals will be to determine the sizes and masses of 
brown dwarf disks. These cannot be measured from exist- 
ing thermal-IR (~3— 10 /um) data, since emission at these 
wavelengths originates from relatively small radii (<1 AU) 
and is optically thick even for minute disk masses. Longer 
wavelength data, namely far-IR (from SIRTF) and sub- 
mm (e.g., from JCMT/SCUBA and SMA) measurements 
are needed. 

7. SUMMARY AND CONCLUDING REMARKS 

We have completed the first systematic survey for cir- 
cumstellar disks around young brown dwarfs and VLM 
stars. We have obtained L'-band imaging for objects spec- 
troscopically classified as having very low temperatures 
and combined these observations with existing JHKs data. 
Our sample comprises most of the published objects with 
spectral types from M6 to M9.5 in IC 348 and Taurus, and 
should be largely free of selection biases for our purposes. 
None appear to be binaries at a resolution of 0'.'4 resolu- 
tion (55-120 AU). Based on current models, our targets 
have masses of ~ 15 to 100 Mjup and ages of <5 Myr. Our 
survey is sensitive enough to detect bare substellar pho- 
tospheres and hence is the first unbiased survey for disks 
around such low mass objects. 

Using the published spectral types and extinctions, we 
determine the intrinsic photospheric Kg — L' (2.0—4.1 /xm) 
colors of our targets and hence if they have IR excesses. 
We find that excesses are very common, occurring in nearly 
80% of the sample. Such excesses are one of the classic sig- 
natures of circumstellar disks around T Tauri stars, and 
there is ample observational evidence and theoretical ex- 



pectation that such disks exists around young stars. Like- 
wise, disks are the natural and most plausible explanation 
for the IR excesses we observe around young brown dwarfs 
and VLM stars. The i'-band emission arises from the in- 
ner regions of the disks, within <0.1 AU, and is sensitive 
to disks with masses of > 10"*" Mq. 

The observed {Ks — L')q excesses are well-correlated 
with Ha line emission, as in the case for T Tauri stars. 
This supports the idea that accretion disks are the com- 
mon origin for both forms of emission. The IR excesses do 
not appear to correlate with any other property of the cen- 
tral sources. The excess frequency and amplitude are in- 
dependent of the mass. However, the coolest objects in the 
sample, spectral type M9, have negligible excesses, indicat- 
ing that their inner disk regions are empty or perhaps the 
objects are diskless. Given the very low estimated masses 
of the M9 objects (< 15 — 20 Mj^p), this is a tantalizing 
finding which should be pursued with a larger sample. 

The T Tauri stars of IC 348 and Taurus exhibit compa- 
rably high disk fractions as the brown dwarfs, and there- 
fore we find that stellar and substellar disks are contem- 
poraneous. Assuming the populations in these regions are 
roughly coeval, brown dwarf disks are at least as long-lived 
as disks around solar-mass stars. The amplitude of the IR 
excesses do not show any trend with age, using the abso- 
lute iC-band magnitude as a model-independent relative 
age scale. Most of the stellar population has an estimated 
age of 1 3 Myr, and assuming coevality, we infer that the 
inner regions of brown dwarf disks do not evolve signifi- 
cantly for the first ^3 Myr. 

Given the low luminosities and weak gravitational po- 
tentials of substellar objects, one expects a priori that 
circumstellar disks around young brown dwarfs will be 
cooler and less luminous than disks around the higher mass 
T Tauri stars. This is confirmed by our measurements: 
the observed {Ks — L')o excesses arc substantially smaller 
than those for T Tauri stars. As a consequence, a con- 
ventional analysis of the disk fraction using only IR colors 
would have missed most of the sources with excesses in our 
sample. While such analyses work well for T Tauri stars, 
in the case of young brown dwarfs the much reduced con- 
trast between the disks and the photospheres is a serious 
impediment. 

We use standard models of fiat, optically thick circum- 
stellar disks to examine the physical properties of the 
brown dwarf disks. The IR excesses can be explained sim- 
ply by reprocessing disks — there is no need for a signifi- 
cant contribution from accretion luminosity. The inferred 
accretion rates are < 10~^ Mq yr~^, or at least an or- 
der of magnitude lower than typical for classical T Tauri 
stars. Nevertheless, the presence of Ha emission strongly 
suggests that at least some accretion is ongoing. The mere 
existence of accreting brown dwarfs at ages of a few Myr 
argues for mass-dependent accretion rates, since brown 
dwarfs with typical T Tauri star accretion rates would not 
remain substellar. 

The maximum observed excesses are less than those ex- 
pected from a disk which extends to the stellar surface — 
this suggests the disks have inner holes. This inference 
is supported by analyzing the observed distribution of IR 
excesses with Monte Carlo simulations of disks viewed at 
random inclinations: models with hole sizes of « 2R^ are a 
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much better match to the observations than models with- 
out holes or those with much larger holes. The presence 
of an inner hole is consistent with a picture where the 
hole arises from the interaction of the inner disk with the 
magnetic field of the central object. Our rough estimates 
based on magnetospheric accretion models find that the 
brown dwarf magnetic fields are on the order of several 
hundred gauss, quite modest compared to measurements 
for T Tauri stars. The inner hole sizes for both T Tauri 
stars and young brown dwarfs are indicative of disk tem- 
peratures hot enough for significant thermal ionization, 
which is needed to permit coupling of the magnetic field 
to the disk. Thus, our inferred hole sizes suggest that the 
magnetic accretion paradigm developed for T Tauri stars 
may extend to substellar masses. This also illustrates one 
way in which brown dwarf disks may be useful laboratories 
for testing our understanding of the physical processes of 
circumstellar disks. 

We have quantitatively examined the possible system- 
atic errors and find that they are not expected to be sig- 
nificant. These include the uncertainties in spectral typ- 
ing and extinction, disk emission contamination, surface 
gravity effects, sample selection bias, and the temperature 
scale for pre-main sequence objects. In addition, most of 
the plausible systematic effects, even if significant, would 
act to increase the disk fraction, which we already find is 
very high. 

Our basic observational findings are that (1) most young 
brown dwarfs have disks, and (2) these disks are contempo- 
raneous with disks around T Tauri stars in the same star- 
forming regions. The latter also demonstrates that brown 
dwarf disks are at least as long-lived as disks around stars, 
assuming that the stars and brown dwarfs are roughly co- 
eval. These observations are naturally accommodated in 
a picture where brown dwarfs are born in a similar man- 
ner as stars — our results offer compelling evidence for a 
common origin for objects from the stellar regime down to 
the substellar and planetary-mass regime. 

Alternative formation scenarios, such as disk-disk colli- 
sions and premature ejection, involve dynamical interac- 
tions in creating brown dwarfs. While specific predictions 
are lacking due to the stochastic nature of these scenarios, 



brown dwarfs formed by collision and/or ejection are gen- 
erally expected to have smaller and less massive disks, and 
consequently shorter disk lifetimes, compared to brown 
dwarfs formed in isolation. This expectation conflicts with 
our finding that brown dwarf disks are at least as long-lived 
as disks around young stars. 

We have found that brown dwarf disks are common. In 
order to better understand the nature of these disks, longer 
wavelength (far-IR, sub-mm, and mm) measurements are 
needed to determine the disk masses and sizes. More mea- 
surements of the accretion rates are sorely needed. The 
lifetimes of disks around brown dwarfs and the disk fre- 
quency at even lower masses (<15 Mjup) are also out- 
standing questions. These studies will be important for 
determining the origin of brown dwarfs and their disks, 
and also for placing these objects in context with our phys- 
ical understanding of the star formation process. Finally, 
the very high disk fraction of young brown dwarfs raises 
the possibility of forming planets around brown dwarfs. 
Such planetary systems would represent a fascinating al- 
ternative to the numerous planetary systems found around 
solar-type stars. 
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APPENDIX 

CALCULATING THE DISK FRACTION IN THE PRESENCE OF OBSERVATIONAL ERRORS 

In the literature, the usual means of computing the disk fraction assumes only Poisson counting statistics contribute to 
the uncertainty in the result. Consider a sample of A'^ objects, each with a color excess A^, which can be derived from a 
single color (as in the case of our work here) or many colors (e.g., a color-color diagram). Typically, one then defines the 
number of sources with disks as 

Ndisk = ^ (Ai > Ao) (Al) 

i 

where the notation represents the number of measurements greater than a specified constant Aq. Many studies choose 
Ao = 0, i.e., any color excess is taken as the sign of a disk (e.g. Lada & Lada 1995; Haisch et al. 2001a). The disk fraction 
and its error are then 

hisk = ^ (A2) 



(^{fdtsk) = ^^^"'^ , (A3) 

respectively. Clearly, this approach neglects the effects of measurement errors: a sample with very large errors in the Ai's 
should lead to a larger a{fdisk) than a sample with small errors. Only in the case of zero measurement errors will (j{fdisk) 
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be given exactly by equation (A3). Hence, the published disk fraction errors are somewhat underestimated. Some studies 
attempt to circumvent the influence of observational errors by setting Aq > 0, typically choosing a value comparable to 
the measurement errors (e.g. Strom et al. 1989; Hillenbrand et al. 1998). However, this approach inevitably excludes some 
objects with genuine disk excesses, and hence the disk fractions are somewhat underestimated. 

We use a simple Monte Carlo approach to account for both Poisson counting errors and measurement errors in deter- 
mining the disk fraction.^ Our observed sample is composed of objects, each with a color index Aj and measurement 
error gj. For each Monte Carlo realization, we generate a set of N objects having color indices 

A^ = Ai + {Gi X e,) (A4) 

where Gi is drawn from a normal distribution of zero mean and unit variance. We then determine the number of sources 

with A^ > 0, denoted as iVjisfe- (The subscript i denotes different objects, and the superscript j denotes different Monte 
Carlo realizations.) For each realization, we assume the resulting number of sources with disks is subject to counting 

statistics: namely, the number of disk sources is described by a Poisson distribution of mean N^^^f., which we denote as 

'^i^disk'- ^disk)- After a large number of realizations (we use 100,000), we sum all the individual Poisson distributions to 
form a probability distribution function for the mimbcr of disk sources: 

PDF(7V,,,,) =Y.V{Nl,^^{Ni^^) . (A5) 
j 

We then compute the disk fraction using the value of N^isk at the peak of the PDF and use the standard deviation of the 
PDF (well approximated by fitting a gaussian) as the error in Ndisk- 

The net result of using this Monte Carlo approach is that our tabulated disk fractions are ^5% smaller than simply 
using the standard approach (i.e., equation (Al) with Aq = 0) and the errors in the disk fraction are f«5% larger. The 
effects are small because of the relatively high S/N of the {Ks~L')o excess measurements (median error of 0.09 mag). 
Note that these percentages should not be used as a general representation for all disk fraction calculations. The result 
of ignoring the measurement errors in Aj will depend on a particular sample's properties, namely the size of the sample, 
the distribution of excesses, and the measurement errors. 
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Fig. 1. — Dereddened absolute /fs-band magnitude and spectral type for our L'-band sample. Magnitude errors are 
comparable to the symbol size, and spectral type errors are ^0.5 subclasses. The main-sequence is shown as a heavy line. 
Dotted lines are model isochrones while solid lines show models of constant mass in units of Mq. The Luhman (1999) 
temperature scale is used to convert the model temperatures to spectral type. Top: Baraffe et al. (1998, 2002) models. 
The isochrones are for ages of 1, 5, and 10 Myr. Middle: D'Antona & Mazzitelli (1997) models. The isochrones are for 
ages of 0.1, 0.5, 1.0, 5.0, and 10.0 Myr. Bottom: Burrows et al. (1997) models. The isochrones are for ages of 0.1, 0.5, 
1.0, 5.0, and 10.0 Myr. 
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Fig. 2. — Top: Observed Ks — L' colors as a function of spectral type. Typical errors in the spectral type are 0.5 
subclasses or less. The heavy line represents the colors of field M dwarfs (see § 2.3 for references). Bottom: Dereddened 
Ks — L' colors. Most of the objects show IR emission in excess of that expected from their photospheres. 
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Fig. 3. — Relation between intrinsic IR excess and equivalent width of Ha emission in Angstroms for the 23 objects 
where such data are available along with Ha upper limits for 2 objects. The two quantities are correlated at the 3a level, 
based on a Spearman rank correlation test. This supports the idea that the optical and near-IR emission both originate 
from the same phenomenon, namely circumstellar accretion disks. 
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Fig. 4. — Histogram of {Ks — L')o excesses as a function of spectral type. The median ±lcr measurement errors 
(±0.09 mag) are comparable to the bin width. The clear bins represent the Taurus sample; the single-hatched bins 
represent the outer IC 348 sample (Luhman 1999); and the double- hatched bins represent the IC 348 core sample 
(Najita et al. 2000). The distribution of excesses among the three earliest spectral type bins (M5.7-M6.4, M6.5-M7.4, 
and M7.5-8.4) are consistent with being drawn from the same parent population based on the K-S test. Hence, there is 
no strong evidence for the IR excesses being dependent on the mass of the central object. 
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Fig. 5. — IR excess as a function of dereddened absolute Xg-band magnitude. At fixed spectral type, Mks is a proxy 
for the relative age, since the evolutionary tracks for young objects are approximately constant in temperature for a 
given mass. No statistically significant correlation with Mks exists for each spectral type bin, suggesting the inner disk 
regions do not evolve substantially over the first ~3 Myr. The symbols represent the different sub-samples and have the 
same meanings as in Figure 1. Errors in M^s are « 0.2 mag, about twice the size of the plotting symbols. 



22 



Liu et al. 





1.2 . 


1 ... 1 . 


■ 1 ■ ■ 


■ 1 ■ ■ ■ 1 ■ ■ ■ .1 ■ ■ ■ 1 




1.0 \ 








X 

1 


CO to 

o o 


/' 

/ 

/ 


□ 






0.4 i 






- 




0.2 i 






SpT = MID ■ 




0.0 / 


/ Ri„= 


1 R 

* 


SdT = ■ 

^ 1 IVI 

SpT = M9 : 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 




1 1 1 1 1 1 


1 L 1 1 1 




1.0 \ 


1 1 1 1 1 1 




1 1 1 1 1 1 1 1 1 .1 1 1 1 1 




0.8 ; 


/ 

/ 






X 

1 


0.6 \ 




□ 




0.4 ; 










0.2 \ 






SpT = Min ■ 




0.0 1 


1 Ri„= 


3R 

* 


SpT = M63 ■ 

SpT = M9 : 

. 1 ... 1 ... 1 ... 1 




1 ... 1 . 


. . 1 . . 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 
Ks-L' 

Fig. 6. — Flat blackbody reprocessing disk models plotted on JHKsL' color-color diagrams. The locus of main sequence 
stars is shown as a solid line and that of giant stars as a dashed line. (See § 4.1 for details.) The two dotted lines 
represent the reddening vector, and the area in between is the locus of reddened stars. Each plot shows models with 
different inner disk radii, and the symbols indicate the spectral type of the central star. The symbols are equally spaced 
in cos 9 from to 1, where 9 is the viewing angle. 
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Fig. 7. — Dereddened Ks — L' colors for our targets (same data and symbols as in Figure 2) compared with face-on 
reprocessing disk models. Disk models with different inner hole radii are labeled in units of the stellar radius, i.e., 
Rin = IR* means there is no hole. The maximum observed IR excesses are less those expected from a disk with no inner 
hole, and suggests Rin ^ 2i?^, is the more likely situation. Notice also that the observed small/non-existent {Ks — L')o 
excesses of the coolest objects, which are spectral type M9, suggest either disks with very large inner holes or no disks 
at all. 
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Fig. 8. - Comparison of tlic observed IR excess distribution with passive disk models possessing different inner radii. 
Disk models with inner radii {Rm) of 1, 2, and 3i?* are shown. For a given disk model, the observed {Ks~L')o excess 
is divided by the maximum possible model excess (i.e., a face-on disk) to form the normalized color excess £{Ks — L')o, 
which is independent of the central object's temperature. The top panels show the observations. Different shadings 
indicate different sub-samples: Taurus {clear bins), IC 348 outer sample {single-hatched bins), and IC 348 core sample 
{double-hatched bins). The bottom panels plot the normalized color excess distribution for a set of disk models seen from 
randomly chosen viewing angles. The observations agree best with the = 2i?^ models, suggesting that disks around 
young brown dwarfs have inner holes. See § 4.2 for details. 
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Fig. 9. — The effect of surface gravity on K — L' colors, based on models by Baraffe et al. (1998, 2002). Dust-free 
("clear") models arc plotted in black, and dusty models in grey. Models with ages of 1 Gyr (higher surface gravity) are 
converted to spectral types using a dwarf temperature scale. Models with ages of 1 Myr (lower surface gravity) use either 
the dwarf temperature scale (upper plot) or the Luhman (1999) pre-main sequence scale {lower plot). For our objects' 
spectral types, the lower surface gravity of young objects is predicted to lead to wO.l mag bluer colors; the observations 
suggest the effect is even smaller. See § 5.3 for details. 
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Fig. 10. — JHKsL' color-color diagram for our sample. The various lines represent the loci of reddened giant and dwarf 
stars (see Figure 6 caption). Median errors are plotted in the lower right. Based on a conventional color-color analysis, 
one would identify only 11 out of 36 objects (those to the right of the dotted reddening band) as having IR excesses. 
However, our analysis incorporating the objects' spectral types shows that in fact many more (31 out of 36) have IR 
excesses: these are plotted as filled circles (•). The majority of sources with disks are missed in the color-color diagram 
because their IR excesses are modest. A similar analysis using only JHKs colors would miss nearly all the objects with 
IR excesses. 
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Table 1 
Infrared Photometry 



Object Spectral Type Ak (mag) L' (mag) E{Ks — L')q (mag) 
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IC 348: outer sample 



L99 205 


M6.0 ± 0.2 


0.13 ± 0.04 


11.43 ± 0.04 


0.38 ± 0.06 


L99 312 


M6.0 ± 0.2 


0.20 ± 0.04 


12.60 ± 0.08 


-0.14 ± 0.09 


L99 325 


M6.0 ± 0.2 


0.19 ± 0.04 


12.27 ± 0.05 


0.39 ± 0.07 


L99 367 


M6.5 ± 0.5 


0.05 ± 0.04 


12.86 ± 0.12 


0.16 ± 0.14 


L99 382 


M6.5 ± 0.5 


0.30 ± 0.04 


12.69 ± 0.07 


0.37 ± 0.10 


L99 407 


M7.0 ± 0.5 


0.40 ± 0.20 


13.77 ± 0.04 


0.11 ± 0.17 


L99 363 


M8.0 ± 0.2 


0.00 ± 0.04 


13.05 ± 0.05 


0.06 ± 0.07 


L99 405 


M8.0 ± 0.2 


0.00 ± 0.04 


13.20 ± 0.07 


0.15 ± 0.09 


Taurus 


V410 Anon 13 


M6 ± 1.0 


0.42 ± 0.06 


10.00 ± 0.04 


0.25 ± 0.11 


V410 Xray 3 


M6 ± 0.2 


0.09 ± 0.06 


9.79 ± 0.04 


0.10 ± 0.07 


MHO 5 


M6.25 ± 0.2 


0.01 ± 0.06 


9.28 ± 0.04 


0.30 ± 0.07 


CFHT-Tau 1 


M7 ± 0.5 


0.30 ± 0.10 


11.14 ± 0.08 


0.02 ± 0.11 


CFHT-Tau 4 


M7 ± 0.5 


0.30 ± 0.10 


9.15 ± 0.03 


0.47 ± 0.09 


CFHT-Tau 2 


M8 ± 0.5 


0.00 ± 0.10 


11.41 ± 0.06 


0.16 ± 0.10 


CFHT-Tau 3 


M9 ± 0.5 


0.00 ± 0.10 


11.63 ± 0.03 


0.00 ± 0.09 



''This object has Ks — L' — 0.23 ± 0.04, much bluer than the other objects in the 
sample. Based on its JHKsL' colors, it appears to be a background late-type giant 
star. It is excluded from the analysis. 
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Table 2 



Fraction of Sources with IR Excesses 



Spectral 


IC 348 


IC 348 


Total 


Type 


core 


core + outer 


(IC 348 + Taurus) 


5.7-6.4 


8/10 


10/13 


13/16 


6.5 - 7.4 


5/6 


8/9 


10/11 


7.5 - 8.4 


4/5 


6/7 


7/8 


8.5-9.4 


1/2 


1/2 


1/3 


all types 


18/23 


25/31 


31/38 




(73% ± 19%) 


(75% ± 17%) 


(77% ± 15%) 



Note. — The percentages of sources with IR excesses and their 
errors were evaluated using a Monte Carlo technique which ac- 
counts for both measurement errors and Poisson counting statis- 
tics. See Appendix A. 



Table 3 

Flat Blackbody Rei'Rocesslxc Disk: MAXLMrii Excess (ln mags) 



Spectral Type A J 


Ail 


Ai^s 


AL' 


AM' 


MO 


0.27 


0.38 


0.63 


1.38 


1 


.78 


Ml 


0.23 


0.35 


0.60 


1.32 


1 


.72 


M2 


0.20 


0.33 


0.56 


1.26 


1 


.66 


M3 


0.18 


0.31 


0.52 


1.20 


1 


.59 


M4 


0.16 


0.28 


0.49 


1.13 


1 


.52 


M5 


0.14 


0.26 


0.45 


1.05 


1 


.44 


M6 


0.12 


0.24 


0.42 


0.98 


1 


.36 


M7 


0.10 


0.22 


0.38 


0.89 


1 


.27 


M8 


0.09 


0.19 


0.34 


0.81 


1 


.17 


M9 


0.08 


0.17 


0.31 


0.73 


1 


.08 



Note. — Transmission profiles for the JHKs bands are 
from 2MASS, while those for the L' and M' bands are 
from the MKO system. 



